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ABSTRACT 



Context. Bipolar outflows from Class protostars often present two components in their CO spectra that have different kinematic 
behaviors: a smooth outflow wing and a discrete, extremely high- velocity (EHV) peak. 

Aims. To better understand the origin of these two outflow components, we have studied and compared their molecular composition. 
Methods. We carried out a molecular survey of the outflows powered by L1448-mm and IRAS 04166+2706, two sources with 
prominent wing and EHV components. For each source, we observed a number of molecular lines towards the brightest outflow 
position and used them to determine column densities for 12 different molecular species. 

Results. The molecular composition of the two outflows is very similar. It presents systematic changes with velocity that we analyze 
by dividing the outflow in three chemical regimes, two of them associated with the wing component and the other the EHV gas. The 
analysis of the two wing regimes shows that species like H2CO and CH3OH favor the low-velocity gas, while SiO and HCN are 
more abundant in the fastest gas. This fastest wing gas presents strong similarities with the composition of the "chemically active" 
LI 157 outflow (whose abundances we re-evaluate in an appendix). We find that the EHV regime is relatively rich in O-bearing species 
compared to the wing regime. The EHV gas is not only detected in CO and SiO (already reported elsewhere), but also in SO, CH3OH, 
and H2CO (newly reported here), with a tentative detection in HCO^. At the same time, the EHV regime is relatively poor in C-bearing 
molecules like CS and HCN, for which we only obtain weak detections or upper limits despite deep integrations. We suggest that this 
difference in composition arises from a lower C/O ratio in the EHV gas. 

Conclusions. The different chemical compositions of the wing and EHV regimes suggest that these two outflow components have 
different physical origins. The wing component is better explained by shocked ambient gas, although none of the existing shock 
models explains all observed features. We hypothesize that the peculiar composition of the EHV gas reflects its origin as a dense wind 
from the protostar or its surrounding disk. 

Key words. Stars: formation - ISM: abundances - ISM: jets and outflows - ISM: individual (IRAS 04166+2706) - ISM: individual 
(LDN L1448) - ISM: molecules - Radio lines: ISM 



1. Introduction 

Outflows from the youngest protostars often display the simplest 
geometry and the highest degree of symmetry, and they are ex- 
pected to reflect more faithfully than other outflows the proper- 
ties o f the still-mysterious driving wind (e.g. jBachiller & Tafallal 
Il999l) . Among the youngest outflows known, one group stands 
out for its pristine appearance in maps and the presence in their 
CO spectra of distinct extremely high velocity (EHV) features. 
These features are secondary emission peaks at the highest out- 
flow velocities that tend to appear symmetrically placed from 
the driving source along the outflow axis. The prototype of this 
group of outflows is the one driven by the 9 L© Class proto- 
star LI 448 - mm, located in the Perse us molecular cloud and first 
identified by iBachiller et al.l (1199 0*). Molecular observations of 
the L1448 outflow have shown that the moving material in this 
system lies in two different components: a pair of conical shells 
and a highly collimated jet. In the spectra, the emission from 
the shells appears in the form of a gradual outflow wing, while 
the emission from the jet arises from the faster EHV feature 
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dBachiUer et al.lll995t lJ0rgensen et alj|2007t iMaurv erani2010l: 
lHiranoetal.1120101) . 

While the shell of relatively slow gas in the L1448 outflow is 
similar to what is seen in more e volved outflows like LI 55 1 (e.g., 
iMoriartv-Schieven et al.l [l987h . and therefore likely composed 
of accelerated ambient gas, the nature of the EHV component 
is still mysterious. The symmetry, high speed, and fragmented 
appearance of the EHV gas initially suggested that it may con- 
sist of a collection of molecular "bullets" ejecte d by the central 
protostar and traveling along the outflow axis (IBachiller et alJ 
1990). Free-traveling bullets with the observed internal veloc- 
ity dispersion of the EHV ga s, however, woul d likely dissi- 
pate quickly along their path (iRicher et al.lll992h . and because 
of this, a number of alternative interpretations to the EHV gas 
have been offered. One possibility is that the EHV gas repre- 
sents shocked ambient material accelerated by an invisible jet, 
as su ggested by observations of SiO emission (Bachilleret al] 
Il99lh and by nume rical simulations of momentum transfer from 
a jet dChernin et alil994) . Alternatively, the EHV features could 
represent internal shocks along a jet owing to variations in 
the velocity of ejection, similar to those invoked to explaiji 
the c hains of HH objects often seen in optical jets (Raga et al 
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Recent high angular resolution observations of very young 
outflows have started to offer new clues to the nature of the 
EHV gas and its relation to the outflow shells. An ideal tar- 
get for this work has been the outflow from IRAS 04166+2706 
(104166 hereafter), a 0.4 L© YSO in Taurus that presents a num- 
ber of similarities with the L1448 outflo w, including a shell- 
plus-jet geometry and bright EHV features (iTafalla et al.ll2004at 
[Santiago-Garcia et al. 2009). High angul ar resolution observa- 
tions of the 104166 outflow by Santiag cvGarcia et al.1 (|2009|) 
show that the shells and the EHV component have very differ- 
ent kinematics. The low- velocity shells seem to arise from am- 
bient gas accelerated by a wide-angle wind, while the EHV gas 
presents a sawtooth velocity field that is inconsistent with an ori- 
gin in shocked ambient material, but is well explained as result- 
ing from a series of internal working surfaces traveling along a 
collimated jet. The 104166 data, therefore, suggests that the out- 
flow consists of two distinct components having different physi- 
cal origin. 

A number of recent theoretical models h ave also suggested 
that outflows consist of multiple components ("Shang et al.'2006t 
[Banerjee & Pudritz 2006; Machidaetal. 2008; Tomidaet^ 
[2010 ). Of particular interest for the study of the EHV gas is the 
so-called "unified model" of .Shang et al., (.2006.) . These authors 
have modeled the interaction between the wind from a proto- 
star and the surrounding dense gas, and suggested that while the 
outflow shells originate from ambient gas that has been shock- 
compressed by a wide-angle wind, the jet-like EHV component 
represents the central and most collimated part of the protostellar 
wind itself. 

Further testing the possible double-nature of the outflow gas 
can be done through chemical analysis. The chemical composi- 
tion of a gas is a sensitive indicator of its thermal history, so if 
the outflow shell and jet-like components have different physical 
origins, this should be reflected in the components having differ- 
ent chemical abundances. The chemical composition of the out- 
flow gas has been studied in a number of sys tems, both towards 
high mass protostars like Orion-lRc 2 (e.g.. iBlake et al.l[l987h 
and towards low mass objects Hk e NGC1333-IRAS 4, LI 157, or 
BHR 71 ( Blake et al. 1995; Bac hiUer & Perez Gutierre3ll997l: 
iGarav et al.lll998l) . These studies have shown that the outflow 
composition is characterized by abundance enhancements of 
species lik e SiO and CH3OH, which a re indicative of shock 
chemistry (Ivan Dishoeck & Blakelll998h . Previous work, how- 
ever, has not included the very young outflows that have EHV 
gas, and as a result, the chemistry of this outflow component has 
not been surveyed in detail and compared with the composition 
of the rest of the outflow. To fill this observational gap, and to 
test the possible double nature of the outflow material, we have 
carried out a systematic molecular survey of the outflow gas to- 
wards the two objects with most prominent EHV components, 
L1448 and 104166. 



2. Target selection and observations 

The intensity of the EHV feature in high density tracers rarely 
exceeds 0.1 K and is often significantly weaker, so a molecu- 
lar survey of the EHV gas requires integration times often in 
excess of 1 hour per transition. This strongly limits the num- 
ber of targets that can be studied during a typical observing 
run, and for this reason, our molecular survey was targeted to- 
wards only two outflow positions, one in L1448 and the other 
in 104166. These positions were selected for having the bright- 
est CO(2-l) EHV component in each outflow, and were de- 
termined during a preliminary search along each outflow axis. 
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Fig. 1. CO(J=2-l) map of the extremely high velocity (EHV) 
component in the 104166 outflow illustrating the location of the 
molecular survey position with a circle of 16" diameter (the typ- 
ical size of the telescope beam width). T he CO map is from the 
Plateau de Bure observations presented in lSantiago-Garcia et al.l 
(2009). 



For the LI 448 outflow, we selected the position with offsets 
(16", -34") with respect to the dri ving source (g( 72 000) = 
3''25"'38^;9, 5(72000) = +30°44^05": TCuriel et al.l[l99(ih . which 
corresponds to the redshifted lobe of the flow. For the 104166 
outflow, we selected the position with offsets (8", 14") with 
respect to the mm-contin uum peak (a( 72000) = 4''1 9"'42'i5, 
5(72000) = +27°13'36"; ISantiago-Garcia et al.ll2009l) . corre- 
sponding to the blue lobe of the outflow. To illustrate its location, 
we present in Fig.[T]the 104166 target position superposed to a 
high resolution CO(2-l) map of the EHV gas 

In addition to the two survey targets, supplementary posi- 
tions of the LI 448 and 104144 outflows were observed for this 
project. After detecting bright SO emission in the EHV compo- 
nents of the L1448 and 104166 outflows, we searched for this 
molecule towards a number of positions along the axis of each 
flow, to asses whether our target positions were unique in pre- 
senting bright SO and to determine the extent of the SO emission 
in the EHV gas. Also, in order to derive beam-dilution factors to 
correct the observed intensities during the multi-line analysis of 
Section 3, we made small CO Nyquist-sampled maps around the 
two target positions. 

All observations presented here were carried out with the 
IRAM 30m radio telescope near Granada (Spain) during ses- 
sions in 2007 June, 2008 August, and 2009 May-June. A num- 
ber of receiver configurations was used to maximize the number 
of molecular lines observed. The list of the main transitions, rest 
frequencies, and related parameters is given in Table[T] For each 
line, two spectrometers were used in parallel: the VESPA auto- 
correlator, configured to provide velocity resolutions of around 
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Fig. 2. CO(J=1-0), CO(2-l), and '^CO(2-l) spectra towards the target positions of the molecular survey in the L1448 (left) and 
104166 (right) outflows. The intensity scale of the spectra has been truncated to better show the two velocity components of the 
outflow, the "wing" and the EHV regime, whose limits are indicated by vertical dotted lines. Negative features near the ambient 
velocity of each source (Vlsr = 5.0 km s ' in L1448 and 6.7 km s ' in 104166) are artifacts of the wobbler switch observing mode, 
and do not affect the outflow emission studied here. 



Table 1. Main transitions observed"'. 



Transition 


Frequency 


A„; 


E„/fc 


g« 




(GHz) 


(s-') 


(K) 




CO(J=1-0) 


115.3 


7.20 10-** 


5.54 


3 


C0(J=2-1) 


230.5 


6.91 10-^ 


16.6 


5 


"C0(J=2-1) 


220.4 


6.08 10-' 


15.9 


5 


SiO(J=2-l) 


86.8 


2.93 10-5 


6.26 


5 


SiO(J=3-2) 


130.3 


1.06 10-* 


12.5 


7 


SiO(J=5^) 


217.1 


5.20 10--* 


31.3 


11 


SiO(J=6-5) 


260.5 


9.12 10-* 


43.8 


13 


SO(Jn=32-2i) 


99.3 


1.13 10-5 


9.23 


7 


SO(Jn=23-12) 


109.3 


1.08 10-5 


21.1 


5 


SO(Jn=43-32) 


138.2 


3.17 10-5 


15.9 


9 


SO(Jn=65-54) 


219.9 


1.34 lO-"* 


35.0 


13 


CH30H(J,=2_i-l_,)E 


96.7 


2.56 10-' 


4.65 


5 


CH,OH(Ji=2o-lo)A 


96.7 


3.41 10-" 


6.97 


5 


CH30H(Jt=2o-lo)E 


96.7 


3.41 10-" 


12.2 


5 


CH30H(Ji=3o-2o)E 


145.1 


1.23 10-5 


19.2 


7 


CH30H(J,=3-i-2_,)E 


145.1 


1.10 10-5 


11.6 


7 


CH,OH(J,=3o-2o)A 


145.1 


1.23 10-5 


13.9 


7 


CH30H(Ja=5-i-4_,)E 


241.8 


5.81 10-5 


32.5 


11 


CH30H(Ji=5o^o)A 


241.8 


6.05 10-5 


348 


11 


HC3N(J=10-9) 


91.0 


5.81 10-5 


24.0 


21 


HC3N(J= 11-10) 


100.1 


7.77 10-5 


28.8 


23 


HC3N(J= 12-11) 


109.2 


1.01 10-^ 


341 


25 


HC3N(J= 16-15) 


145.6 


2.42 lO-t 


59.4 


33 


CS(J=2-1) 


98.0 


1.68 10-5 


7.06 


5 


CS(J=3-2) 


146.7 


6.07 10-5 


141 


7 


CS(J=5^) 


244.9 


2.98 10-^ 


35.3 


11 


HCO+(J=1~0) 


89.2 


419 10-5 


4.28 


3 


HCN(J=1-0) 


88.6 


2.41 10-5 


4.26 


3 


H2CO(J;f„;,,=2o2-loi) 


145.6 


7.81 10-5 


10.5 


5 


S02(J;f„;,,=3i3-2o2) 


104.0 


1.01 10-5 


7.75 


7 


H2S(Js:„a:^ = 1io-1oi) 


168.8 


2.68 10-5 


27.9 


9 


SiS(J=5^) 


90.8 


1.19 10-5 


13.1 


11 


HNC(J=1-0) 


90.7 


2.69 10-5 


4.35 


3 


(1) Molecular parameters from the Cologne Database for Molecul 


Spectroscopy (http://v™w. astro. 


xmi-koeln.de/cdms/), s 


iMulleretal. (2001, 2005^ 











1 km s-', and either the 1 MHz filterbank or the wider (4 MHz) 
WILMA autocorrelator to provide extra frequency coverage. 
Flat baselines were achieved using wobbler-switching mode, 
chopping with a frequency of 0.5 Hz between the source position 
and a reference position 220" away in azimuth. The relatively 
small off'set between source and reference position caused some 
contamination in the spectra at the ambient cloud velocity, but it 
did not affect the outflow emission studied here. As a result of 
the flat baselines, only a zeroth-order polynomial was subtracted 
from the spectra in the off-line data reduction. 

During the observations, the pointing was checked and cor- 
rected by making cross scans on nearby continuum sources ap- 
proximately every two hours, and it was usually below 3" rms. A 
number of larger pointing excursions occurred during the 2008 
run, and observations with suspicious pointing errors were dis- 
carded. Calibration was achieved by observing every 10 minutes 
a sequence of blank sky, an absorber at room temperature, and 
a cold load, and the calibrated T^ scale from the telescope was 
converted into the r^b scale using the standard beam efficiency 
factors provided by the telescope system. Further reduction and 
analysis of the data was carried out using the GILDAS software 
dhttp : //www ■ iram. fr/IRAMFR/GILDASl . 

3. Molecular survey results 

3. 1 . CO data: wing and EHV regimes 

Figure |2] shows the CO(l-O), CO(2-l) and '^CO(2-l) spectra 
towards the two target positions of our survey. As can be seen, 
each CO spectrum presents a broad wing component that starts 
at the ambient cloud velocity (Vlsr = 5.0 km s-' in L1448 and 
6.7 km S-' in 104166) and extends towards the red in the case of 
L1448 and towards the blue in 104166. This wing component de- 
creases gradually in intensity without any sign of discontinuity 
until it merges at high velocities with the bright EHV compo- 
nent. This EHV component appears in the spectrum as a sep- 
arate secondary peak of intensity 0.5-1 K and centered around 
Vlsr = 70 km s-' in L1448 and -35 km s-' in 104166. In L1448, 
the EHV component is detected just above the noise level in the 
'^CO(2-l) spectrum, while in 104166, the EHV feature can only 
be seen in the main isotopologues. 
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To carry out our abundance analysis, we need to define ve- 
locity limits for the wing and EHV components of each outflow, 
and we use for this the CO spectra in Fig.|2] For the EHV regime, 
we define its limits from the velocities at which this feature dom- 
inates the emission in the spectrum: Vlsr from 55 to 85 km s"' 
in L1448 and -26 to -43 km s ' in 104166 (vertical dotted fines 
in Fig. |2]l. For the wing regime, we take as high velocity limit 
the lower value of the EHV component, as there is no clear gap 
between the two emissions in the spectra. Selecting the low ve- 
locity limit for the wing regime requires an additional consider- 
ation, as the wing emission merges smoothly with the optically 
thick ambient component. To guarantee that our wing contribu- 
tion is optically thin in CO, we have calculated the ratio between 
the CO(2-l) and '^CO(2-l) spectra and determined the velocity 
at which the ratio becomes higher than 50. This ratio is approx- 
imately equal to the expecte d '^C/'^C ab undance ratio in the lo- 
cal ISM (k 57. ILanger & Penziasl|19 90), and is comparable to 
the maximum value measured with our data due to the limited 
S/N ratio of the '^CO(2-l ) spectrum. According to our estimate, 
it is reached at Vlsr - 10 km s"' in L1448 and at 4 km s"' in 
104166, which from now on will be used as lower velocity limits 
for the (optically thin) wing regime. 

Although the above velocity limits help separate the outflow 
into its two main components, velocity is not the only prop- 
erty that defines the wing or the EHV regime. Interferometer 
maps of the CO emission show that each regime corresponds 
to a different part of the outflow, with the wing emission aris- 
ing from a shell an d the EHV component lying along a jet (e.g. 
ISantiaeo-Garcia e t al. 2009). These two components may there- 
fore appear to overlap in velocity towards certain directions, but 
this does not necessarily imply that they spatially coexist in the 
outflow. Indeed, we will see below that the wing component 
seems to continue to higher velocities inside the EHV regime. 
Such an extension is very weak, and can be usually ignored when 
estimating the EHV contribution. It can however be detected in 
species with weak EHV emission, like HCN and possibly CS, 
and in those cases, it needs to be handled with some care. 

3.2. Overview of the molecular survey 

Fig |3] shows the spectra of the different high density tracers de- 
tected towards the L1448 and 104166 outflow target positions. 
As can be seen, a number of molecules have detectable emission 
in one or both outflow regimes, and overall, the spectra of most 
species resemble the CO spectra discussed in the previous sec- 
tion. Bright emission in the wing regime is seen, among others, 
in the lines of HCN, SO, and in the low J transitions of SiO, 
while weaker but still noticeable wing emission can be seen in 
the spectra of CH3OH, CS, and HCO+. 

As Fig[3]also shows, several species present bright emission 
in the EHV regime. The SiO lines display the brightest EHV fea- 
tures, and have peak intensities of ~ 100 mK. EHV SiO emis- 
sion has previously been detected both in L1448 and 104166, 
and this molecule has been so far the only known tracer of 
flie EHV gas apart from CO (Bachill er et al.lll99ll: fNisini et al.1 
l2007t ISantiago-Garcia et al.l |200 9). As our observations indi- 
cate, SiO is not the only dense gas tracer detectable in this out- 
flow component. The SO molecule, for example, also presents 
significant EHV emission in both L1448 and 104166. EHV com- 
ponents were detected in SO(Jn=32-2i) and SO(Jn=43-32) with 
an intensity of x 50 mK, or about half of that of SiO, and in 
the brighter L1448 outflow, the higher excitation SO(Jn=65- 
54) transition was also detected at a 4cr level (not shown). 
Additional species with definitive detection of EHV emission 



are CH3OH (L1448 and likely 104166) and H2CO (LI448 only). 
They present intensities on the order of 10 mK, and their discus- 
sion is deferred to section 3.4, where a number of weak detec- 
tions and upper limits for the EHV gas will be further studied. 

Before starting our analysis of the outflow chemical compo- 
sition, we note that two main results from our survey can already 
be inferred from the direct inspection of the spectra in Fig . [3] The 
first result is that despite their factor of 10 difference in energet- 
ics, the L1448 and 104166 outflows present very similar excita- 
tion and chemical composition in both the wing and the EHV 
components. This can be noticed from the similar trends seen in 
Fig.[3]towards the two outflows. For example, in SiO, both out- 
flows present an EHV emission with all transitions up to J=5^ 
having similar intensities (» 100 mK) while, again in both out- 
flows, the wing emission fades quickly in the high-J lines. Also, 
both outflows present a combination of relatively weak CS and 
HCO^ wings and very prominent HCN wing. This all suggests 
that the pattern of abundances in the two objects is very similar 

A second result that can be derived from an inspection of 
the line profiles is that the wing and the EHV components dif- 
fer significantly in their chemical composition. This can be in- 
ferred from a lack of correlation between the intensity of the 
wing emission and the detection or not of an EHV feature. An 
extreme example of this behavior is the HCN(l-O) line, which 
has stronger wing emission than transitions of similar Einstein-A 
coefficient, like those of SiO and SO, but that has no detectable 
EHV secondary feature in the spectrum (the HCN emission in 
this range seems only the continuation of the wing, see section 
3.4). Such a different behavior of the two outflow components 
is a first indication that they have undergone different chemical 
processing, which is a hint of a possible different physical origin 
of the wing and the EHV gas. Before further discussion, how- 
ever, we need to investigate the internal chemical structure of 
the wing component. 

3.3. Two chemical regimes in the outflow wing component 

The outflow wing appears in the CO spectra as a single compo- 
nent with no internal structure, but in species like SiO, CS, or 
H2CO, the wing presents a more complex shape. To illustrate 
this behavior, we present in Fig. |4] spectra for a few representa- 
tive species together with scaled-down versions of the CO line 
profile. As can be seen, especially in the higher S/N L1448 data, 
the lines present different wing shapes, and only the SO spec- 
trum resembles the CO profile over the whole wing regime. The 
CS and SiO profiles, present two different wing sections: a low 
velocity part in which the intensity lies below the (scaled-down) 
CO spectrum and a high velocity part in which the intensity fol- 
lows the CO line shape (a vertical dashed line has been drawn to 
separate these two wing sections). In contrast with SO and SiO, 
H2CO presents a different (and unique) behavior Its line profile 
seems to match the CO spectrum in the low velocity section of 
the wing while it drops below CO at higher velocities. As Fig.|4] 
shows, despite the variety of line profiles and wing sections, the 
behavior of each species is remarkably similar in both outflows. 
This suggests that the line changes we see in the figure are more 
related to the peculiar chemistry of each species than to the de- 
tails of the outflow. 

Another remarkable property of the line profiles in Fig.|4]is 
the transition between the two wing sections. For a given out- 
flow, this transition occurs at the same velocity for all species, 
and often in a relatively sudden way. In L1448, both CS and 
SiO spectra have a » 5 km s^'-wide velocity section just be- 
low Vlsr - 20 km s"' where the wing shape reverses, and the 
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Fig. 3. Spectra from the high density tracers with detected outflow emission in the L1448 (left) and 104166 (right) survey targets. 
As in Fig. |2] the vertical dotted lines mark the boundaries of the wing and EHV outflow regimes. Note the prominent EHV features 
in the SiO and SO lines, and the prominent wings of lines from other species like HCN. Note also how the outflow regimes partly 
overlap in the multicomponent spectra of CH3OH. The analysis used to disentangle these components is described in Appendix C. 



intensity of the emission increases slightly with velocity before 
reaching the faster wing section where it matches the CO pro- 
file. It is unlikely that this dip in the wing emission arises from 
self-absorption by ambient gas because the dip is red-shifted by 
10 km s"' with respect to the ambient velocity. The dip is also 
unlikely to arise from contamination by the reference position 
because it was seen under different observing conditions and in 
different epochs. It must represent a real feature of the emission. 
In 104166, no sharp transition is seen between the two wing sec- 



tions, but the CS and SiO spectra seem to match the fast part of 
the CO wing at the same velocity near Vlsr = -7 km s"'. 

The above differences in the shape of the profiles indicate 
that for most molecules, the ratio of intensity with respect to 
CO changes with velocity over the wing regime. Barring exci- 
tation and optical depth effects (discussed in the following sec- 
tions), such a change in the intensity ratio must correspond to 
a change in the ratio of the abundance with respect to CO. The 
wing regime, therefore, seems not made of gas having uniform 
composition, but to consist of at least two regimes with different 
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Fig. 4. Comparison of spectra from selected dense gas tracers (black) with scaled-down versions of the CO(2-l) spectrum (red) 
to illustrate the existence of two different chemical components in the wing regime of the L1448 and 104166 outflows. Note how 
in the SO, CS, and SiO spectra, the fast wing (f-wing) component matches well the CO spectrum, while the slow wing (s-wing) 
component follows a different trend: in SO, it matches the CO spectrum, while in CS and SiO, it lies below it. In contrast, the H2CO 
wing matches the CO s-wing but lies below the f-wing at higher velocities. This different behavior of the spectra in the s-wing 
and f-wing regimes suggests that the composition of these two parts of the outflow follow different trends and need to be studied 
separately 



Table 2. Velocity ranges of the slow wing (s-wing), fast wing 
(f-wing), and extremely high velocity (EHV) regimes. 



Target position 



s-wing 
(kms-') 



f-wing EHV 

(kms"') (kms"') 



L1448 (16", -34")"> 
104166 (8", 14")<2) 



[10, 20] 
[-7, 4] 



[20, 55] 
[-26, -7] 



[55, 85] 
[-43, -26] 



(1) Offsets with respect to Q'(y2000) = 3''25'"38^;9, (5(72000) = 
+30° 44' 05"; (2) Offsets with respect q'(72000) = 4'T9'"42!5, 
(5(72000) = +27° 13' 36". 



chemical properties. To include this behavior in our abundance 
analysis, from now on we will consider the two wing regimes 
separately, and we will refer to them as the slow wing or s-wing 
and the fast wing or f-wing. Their velocity limits, and those of 
the EHV component, are indicated by vertical dotted lines in 
Fig. |4] and summarized in Table |2] 



3.4. Molecular emission in the EHV regime 

The emission of molecules in the EHV regime deserves spe- 
cial attention. Before our survey was carried out, only CO and 
SiO had been detected in this outflow component, but as shown 
in section 3.2, SO also presents bright EHV components in at 
least two transitions towards both L1448 and 104166. The SO 
molecule is therefore a relatively bright tracer of this fastest out- 
flow component, and in Appendix A we present the results of a 
small SO survey along the axes of the L1448 and 104166 out- 
flows illustrating the widespread nature of the EHV SO emis- 
sion. 



In addition to CO, SiO, and SO, our observations revealed 
weak (~ 10 mK) high-velocity emission in a number of addi- 
tional species, some of them presenting a distinct EHV feature 
while others of a less clear nature. This is illustrated in Fig. |5] 
which shows an expanded view of the spectra for all species 
with detectable emission at high velocities, together with scaled- 
down versions of the SO(Jn=32-2i) line to indicate the expected 
position of the EHV regime. (For the CH3OH and HCN lines, the 
SO models are in fact templates of the multi component spec- 
trum and were created using the method described in Appendix 
C.) Given the weak emission of the features, the variety of line 
shapes, and the possible contamination from lines of additional 
molecules, we will consider each line separately. 

As Fig. Sshows, both the CH30H(2i-li) and CH30H(3i- 
2k) spectra towards the L1448 outflow present distinct secondary 
peaks in the velocity range expected for the EHV component 
(Vlsr - 60-80 kms"'.) These peaks match closely the EHV fea- 
tures predicted by the SO template, indicated in the figure by a 
red line plotted below the CH3OH spectrum. For CH30H(3t- 
2k), the SO template predicts a significantly broader EHV com- 
ponent due to the overlap of the contributions from the E and A 
CH3OH species, and indeed the observed CH30H(3i-2i) EHV 
feature is significantly broader than in CH30H(2i-li). We con- 
sider this double match between observations and SO templates 
a definitive detection of CH3OH in the EHV gas of the L1448 
outflow. 

The situation in 104166 is harder to evaluate due to a signif- 
icantly (factor of 2) weaker CH3OH signal and to the expected 
overlap in CH30H(3i-2,i) between the EHV feature and a group 
of weak (mostly k - 2) CH3OH ambient-velocity components 
whose position is indicated in the figure by a horizontal bar. As 
Fig. |5] shows, both CH3OH spectra present significant emission 
at the expected velocity of the EHV gas, but the lack of a distinct 
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Fig. 5. Evidence for EHV emission in CH3OH and H2CO towards the L1448 outflow, together with possible detections in HCO^, 
CS, and HCN. The black histograms represent the data and the red histograms are scaled-down SO(32-2i) spectra to illustrate the 
expected location of the EHV peak (near Vlsr = 70 km s ' in L1448 and near -35 km s ' in 104166). For CH3OH and HCN, the 
red histogram is a synthetic spectrum derived using the SO(32-2i) line as a template (see Appendix C for further details). Note 
the weak EHV peaks and the gradual transition from well defined secondary components in CH3OH and H2CO to more wing-like 
features in HCO^, CS, and HCN. Several unrelated lines are indicated, and the horizontal lines in the CHj,OH(3i!-2ii) spectra mark 
the position of weak k - 2 components that may contaminate the EHV feature in 104166. See text for a discussion of the reliability 
of each detection. SO was chosen as a template due to its similar behavior to CO and higher sensitivity to dense gas. 



secondary peak isolated from the rest of the outflow wing makes 
this detection less clear than that in the L1448 outflow. 

Another convincing detection of an EHV feature can be seen 
in the H2CO spectrum towards the L1448 outflow. Fig. |5]shows 
that the H2CO(2o2-l{)i) line presents a well-defined secondary 
peak at the expected velocity for the EHV gas as indicated 
by the SO profile, again shown in red below the H2CO spec- 
trum. This H2CO feature stands out in the spectrum as a distinct 
component, and is well separated from the brighter peak near 
Vlsr > 100 km s"' that corresponds to the ambient and outflow 
wing contributions from HC3N(16-15). We thus consider this 
feature as a clear detection of H2CO in the EHV component of 
the LI 448 outflow. 

As with CH3OH, no clear H2CO EHV feature can be seen 
in the weaker 104166 spectrum. Our data only provides an up- 
per limit to the H2CO emission in this outflow component, and 
a deeper integration is needed to test whether there is any differ- 
ence between the presence of H2CO in the EHV component of 
the two outflows. 

The third species with a likely EHV feature in the spectrum 
is HCO+. As Fig.|5]shows, the HCO+(1-0) fine from the L1448 
outflow presents an almost continuous wing with a hint of a rela- 
tive maximum near the EHV velocity range. In the 104166 spec- 
trum, on the other hand, HCO^(l-O) presents a distinct emission 



peak at the velocity expected for the EHV regime, suggesting 
that at least part of the HCO^ emission at high velocities truly 
arises from the EHV gas. These two HCO^ spectra illustrate a 
trend, further seen in CS and HCN, of increasing confusion be- 
tween the EHV feature and what seems to be a high-velocity 
continuation of the wing emission. As mentioned in Sect. 3.1, 
the wing emission in an outflow probably does not end at the ve- 
locity at which the EHV component starts to dominate the spec- 
trum, but it continues to higher velocities overlapping with the 
EHV gas. In the lines of CO (also SiO and SO), the EHV emis- 
sion is brighter than the wing, so it dominates the spectrum. In 
HCO+, CS, and HCN, the weak EHV emission barely (if at all) 
stands out over the wing emission, and its separation from the 
fastest wing is ambiguous. 

The last two species that we consider are CS and HCN. Their 
spectra present statistically significant emission at velocities of 
the EHV regime (bottom two panels in Fig. [Sj, but the nature 
of this emission is unclear. The CS(2-1) spectrum in both the 
L1448 and 104166 outflows presents rather flat emission in the 
expected range of the EHV component, with a possible hint of 
a secondary peak in the L1448 outflow. The 104166 spectrum 
may be slightly contaminated by a pair of C3H lines, so its shape 
should not be used to judge a possible EHV peak. The HCN(1- 
0) spectrum, on the other hand, presents smooth wing-like emis- 
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sion at the velocities expected for the EHV component (as pre- 
dicted by the SO template), suggesting that in this molecule the 
wing emission overwhelms any possible contribution from the 
EHV gas. Thus, while it is possible that the EHV component 
emits non-negligibly in CS(2-1), we consider that most of the 
high velocity HCN(l-O) emission arises from the extended wing 
component and not from the EHV gas. 

In summary, our observations show definite emission from 
the EHV component in CH3OH and H2CO towai'ds the L1448 
outflow, likely emission in HCO^ for both outflows, possible 
emission in CS, and unlikely emission in HCN. The results to- 
wards 104166 are compatible with those towards L1448. Given 
the several hours of integration needed to obtain the spectra of 
Fig-El it seems unlikely that a more definitive conclusion for the 
marginal cases can be reached in the foreseeable future. 

3.5. Column density estimates 

We start our abundance analysis by estimating the column den- 
sity of each species in our survey. As a first step, we integrate the 
observed intensity inside each of the three velocity ranges into 
which we have divided the outflow emission (s-wing, f-wing, 
and EHV, see Table |2]i. For species like CO, SiO, SO, CS, and 
CH3OH, we have observed more than one transition, so we can 
use a population diagram analysis to check whether the molecu- 
lar energy levels are in LTE, and if so, t o derive an excitation 
temp erature and a total column density (JGoldsmith & Langed 
11999 ). Throughout this analysis, we assume that the outflow 
emission is optically thin. For the CO lines the condition is guar- 
anteed by the comparison with the '^CO data (section 3.1), and 
for other species, the data already provides a strong hint of the 
emission being thin because the observed intensities are on the 
order of 0. 1 K, while the excitation temperatures are close to 7 K 
(see below). 

Before converting the observed line intensities into col- 
umn densities, a number of corrections are needed. A first cor- 
rection is required to compensate for the fact that transitions 
with different frequencies were observed using different beam 
sizes. As interferome ter maps show (e.g., Dutrev et al. 1997; 
ISantiago-Garcia et al.l 12009). the outflow emission is generally 
compact, especially in the EHV regime, so the smaller beam size 
of the high frequency transitions is expected to enhance their 
emission. As these lines usually correspond to high excitation 
levels, this effect can introduce a systematic bias in the column 
density estimate, and needs to be corrected. For this, we have 
used high-resolution CO data to model the dependence of the ob- 
served intensity with beam size in each of the outflow regimes, 
and we have derived correction factors to convert the emission 
of each transition into the expected value for a 16" beam (the 
approximate mean value of our beam sizes). The good match 
between line shapes of the high density tracers and CO (Fig.|4]i 
suggests that these CO-based correction factors can be applied 
to the different species in our survey. A full description of the 
procedure is presented in Appendix B. 

A second correction is required to disentangle the overlap 
between multiple lines in the spectra of HCN and CH3OH. In 
these spectra, the ambient lines are usually well separated, but 
the outflow regimes are not. Fortunately, the overlap between 
outflow components is only partial, and the outflow contribu- 
tion to each individual line seems to always produce the same 
type of profile. Thus, we have modeled each multi-component 
spectrum as consisting of partly overlapping lines all having the 
same spectral shape and each centered at the velocity given by 
their rest frequency. As templates for these components, we have 
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Fig. 6. Population diagrams for CO, SiO, and SO in L1448 and 
104166 illustrating the results from the excitation analysis. Red 
triangles, blue squares, and black circles represent s-wing, f- 
wing, and EHV values, respectively. The lines show the results 
from the least-squares fit described in the text and summarized 
in Table [3] The SO data have been multiplied by 50 for presen- 
tation purposes only. 



used observed lines, usually from SO and SiO, and we have de- 
rived their relative intensities from a fit to the full spectrum (see 
Appendix C for further details). 

A final consideration in the column density analysis is 
needed because the excitation temperatures derived for the out- 
flow gas are generally low due to subthermal excitation (ss 7 K), 
so that the common assumption that T^^ is much higher than Thg, 
the cosmic background temperature, is not guaranteed. In this 
situation (and assuming that the line is thin), the column density 
of the upper level of one transition (A^„) is related to the observed 
integrated intensity W by 



N, 



hc^A„i 



Jiv/kT„ 



1 



n-1 



i,hvlkT,,g 



1 



where k is Boltzmann's constant, v the line frequency, h Planck's 
constant, c the speed of light, and A„/ the Einstein A-coefficient. 
The term in square brackets, ignored in the standard (high- 
temperature) population diagram analysis, makes the column 
density estimate depend on T^^, which is not known until the 
analysis has been performed. Fortunately, this term is usually 
close enough to 1 that A^„ can be solved iteratively. In our com- 
putation, we first assume that the term is equal to 1 and derive 
Tex using the standard population analysis. Then, we iterate the 
LTE estimate recalculating A^„ with the r^ , value from the pre- 
vious iteration. Convergence to a level of lO"'* in r^ , is usually 
reached in 4-5 iterations. 

As a result of our LTE analysis, we have derived T^x and 
column density {Niot) values for each outflow component in all 
species with more than one observed transition. These values are 
presented Table[3] and population diagrams for CO, SO, and SiO 
are shown in Fig.|6]for illustration. As the table and figure show, 
the Tex values for CO are close to 15-20 K in both outflows. 
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Table 3. Molecular column densities and excitation temperatures derived towards the L1448 and 104166 outflow targets 

L1448 



MOLEC 



s-wing 
Nj- (cm--) 



T„ (K) 



f-wing 



T„ (K) 



EHV 



Nr (cm-^) 



T.. (K) 



CO 


1.0 10'" ± 1 10"* 


18.6 ±0.3 


7.4 10'^ ±2 lO'"* 


15.6 ±0.5 


5.7 10'^ ±2 lO''' 


11.1 ±0.3 


SiO 


2.4 10'- ±2 10" 


5.6 ± 0.3 


4.7 10'2±2 10" 


7.1 ±0.3 


6.4 10'2±2 10" 


7.0 ±0.3 


SO 


1.8 10'3±6 10" 


7.2 ±0.1 


1.4 10'3±9 10" 


7.9 ±0.3 


1.1 io'-''±i 10'2 


8.3 ± 0.5 


CS 


5.5 10'- ±2 10" 


6.3 ±0.1 


8.5 10'2±3 10" 


6.1 ±0.2 


7 10" ±3 10" 


7 ± IC) 


CH^OH<2) 


7.9 10" ± 1.1 10'^ 


10 ±1 


5.1 10'3±1.5 10'^^ 


11±4 


1.6 10'3±5 10'2 


10 ±3 


HC3N 


2.0 10'- ±5 10" 


14 ±2 


2.4 10'^ ± 1.4 10'- 


11±3 


< 1 10'2 


7± 1<" 


HCO+ 


9.3 10" ±2 10'" 


7±1"' 


7.0 10" ±3 10'" 


7±1<" 


5.0 10" ±4 10'" 


7±1"' 


p-H,CO 


1.9 10'- ±3 10" 


7 ±10 


7.4 10" ±1 10" 


7±1<'> 


2.3 10" ±3 10'" 


7±1<" 


HCN 


9.7 10'- ± 1 10'2 


7±1(') 


1.6 10'3±3 10" 


7±1<" 


1.2 10'2±2 10" 


7±1<" 


SO. 


3.4 10'2±6 10" 


7±1"' 


< 2 10'2 


7±1<" 


< 3 10'2 


7± 1<" 


SiS 


< 5 10" 


7 ±10 


< 2 10'^ 


7±1<" 


< 1 10'2 


7± 1(" 


H,S 


< 3 10'^ 


7±1<" 


<5 10'- 


7±1<" 


<4 10'2 


7± 1"> 



104166 



MOLEC 



s-wing 
Nj- (cm-2) 



T„ (K) 



f-wing 



T„ (K) 



EHV 



Nr (cm-2) 



T„ (K) 



CO 


6.5 10'^ ± 1 10'-* 


24.5 ± 0.6 


1.5 10'^ ± 1 10'^ 


17 ±1 


3.2 10'^ ± 1 10'^ 


18±1 


SiO 


1.2 10'- ±2 10" 


7.3 ±0.7 


1.1 10'2±1 10" 


8.0 ±0.6 


4.0 10'2±1 10" 


9.6 ±0.1 


SO 


1.2 10'^^ ±8 10" 


6.0 ±0.2 


3.6 10'2±2 10'- 


5.0 ±0.8 


7.3 10'2±1 10'^ 


6.6 ±0.4 


CS 


2.1 10'- ± 1 10" 


5.6 ± 0.2 


1.4 10'2±3 10" 


4.6 ± 0.5 


5.4 10" ±6 10'" 


7±1<" 


CH,OH <2) 


5.2 10'^^ ± 8 10'2 


11±1 


8.1 1012 ±2 10'2 


8±2 


6.3 10'2±1.8 10'2 


7±2 


HC3N 


< 2 10'- 


7±1<" 


<5 10" 


7±1<" 


<5 10" 


7± 1<'> 


HCO+ 


2.6 10" ±2 10'" 


7±1(') 


1.9 10" ±3 10'" 


7±1<" 


3.0 10" ±2 10'" 


7± 1<" 


p-HoCO 


1.8 10'- ±2 10" 


7±1(" 


< 6 10'" 


7±1<" 


<5 10" 


7± 10 


HCN(^^' 


3.2 10'- ±3 10" 


7 ± 1'" 


3.4 10'2±3 10" 


7±1<" 


2.8 10" ±3 10'" 


7±1"' 


SO9 


< 2 10'2 


7±1(') 


< 1 10'- 


7±1<" 


< 1 10'2 


7± 1<" 


SiS 


< 5 10" 


7 ± 1<" 


< 1 10'2 


7±1<" 


<7 10" 


7± 10 


H2S 


< 5 10'- 


7±1<" 


<2 10'^ 


7±1<" 


<3 10'^ 


7± 1'" 



Notes: (1) r„ = 7 ± 1 K assumed; (2) The E and A-forms of CH3OH were analyzed separately and the column density in the table is the sum 
of the two results (r^j is the mean value). The E/A ratio was found » 1.2 in both sources and for all outflow regimes; (3) 10% uncertainty 
assumed for 104166 estimate of HCN column density. 

Table 4. CO-normalized abundances (xlO^) and enhancement factors with respect to dense core abundances''^ 



L1448 



104166 



MOLEC 



s-wing 

vCO f 



f-wing 

yCO f 



EHV 



X^ 



s-wmg 

vCO f 



f-wing 

vCO f 



EHV 



XV 



SiO<2) 


24 


4.0 10' 


64 


1.1 10^ 


11 


1.9 10^ 


1.8 


3.1 10' 


7.3 


1.2 10-* 


12.5 


2.1 10* 


SO 


18 


2.7 10^ 


19 


2.9 10^ 


19 


2.9 10- 


18 


2.8 10^ 


24 


3.7 10^ 


23 


3.5 10^ 


CS 


5.5 


7.9 


11 


17 


1.2 


1.7 


3.2 


4.7 


9.3 


13 


1.7 


2.4 


CH30H 


79 


5.7 10^ 


69 


5.0 10^ 


28 


2.0 10^ 


80 


5.8 10^ 


54 


3.9 10^ 


20 


1.4 10^ 


HC3N 


2.0 


13 


3.2 


21 


< 1.8 


< 11 


<3.1 


<20 


<3.3 


<22 


< 1.7 


< 10 


HCO+ 


0.93 


2.0 


0.95 


1.9 


0.88 


1.8 


0.4 


0.8 


1.3 


2.5 


0.9 


1.9 


p-H2CO<3' 


1.9 


27 


1.0 


14 


0.40 


5.5 


2.8 


39 


<0.4 


<5.7 


< 1.4 


<22 


HCN 


9.7 


9.1 


22 


20 


2.1 


2.0 


4.9 


4.6 


23 


21 


0.88 


0.82 



Not es: (1) X^° = lO^'x Nt/N(CO), and /„,/, as defined in Section 4.3; (2) The assumed SiO abundance in cores is the upper limit suggested by 
IZiurvs et alj ( ll98^ : 5 10"'^; (3) An ortho-to-para ratio of 3 was assumed for H2CO. 



while the Tg^- values for the high-density tracers are lower than 
10 K and usually close to 7-8 K, also in both outflows and in 
all the velocity regimes. This difference in T^v between CO and 
the other tracers indicates that the emission we observe in all the 
outflow regimes arises from gas where the high-density tracers 
are subthermally excited while CO is probably thermalized (if 
CO were subthermal, then T^ i for the high density tracers would 



be close to the cosmic background temperature). Such excitation 
conditions are satisfied if the emitting gas has a kinetic temper- 
ature of around 20 K and a density close to 10^ cm"^, condi- 
tions which mat ch closely those p redicted by realistic models of 
post-shock gas ("Bergin et al.lll998 ). A more detailed analysis of 
the excitation conditions in the outflow gas for a larger sample 
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Fig. 7. Ratio between the CO-normalized abundances derived 
for L1448 and 104166 in each outflow regime (symbols as in 
Fig. |6]l. Note how the points cluster around a value of 1. Only 
detections have been used. No HC3N data are presented due to 
its non detection in the 104166 outflow (at levels consistent with 
an abundance ratio of 1, see Table|3]l. 
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Fig. 9. Comparison between the HCN(l-O) spectrum in L1448 
and a synthetic spectrum generated using SiO(2-l) as a tem- 
plate. Note the good match between the two histograms in the f- 
wing and the large difference in the EHV regime. This change of 
behavior indicates a higher (factor of 20) change in the SiO/HCN 
abundance ratio between the f-wing and the EHV regime. 



of objects will be presented in Santiago-Garcia et al. (2010, in 
preparation). 

The similar Tg^ values shown in Table |3]for the high-density 
tracers in the outflow suggest that for those species for which 
only one transition was observed in our survey (e.g., HCO^, 
H2CO, HCN), we can safely assume a single excitation value in 
the column density calculation. Based on the SiO, SO, and CS 
estimates, we take as default value Te^ = 7 + 1 K, and use it to 
derive column densities for the species with only one observed 
transition. For species with no detection, we use the spectrum 
noise level together with the default Tg^ value to derive upper 
limits to their column density. Table|3]summarizes all the results. 

4. Analysis of the outflow abundances 

4.1. The similar composition of the L1448 and 104166 
outflows 

As Table |3] shows, the CO column densities in the LI 448 out- 
flow are significantly higher than in 104166. The ratio between 
the two is 1.5 in the s-wing regime, 4.9 in the f-wing regime, 
and 1.8 in the EHV regime, and these numbers are in general 
agreement with the fact that the L1448 outflow is more massive 
and energetic than that of 104166. To compensate for this dif- 
ference, and to properly compare the abundance of the different 
molecular species in the two sources, we normalize each column 
density by the CO value in the same regime, working from now 
on with abundances with respect to CO. Such CO-normalized 
abundances (referred to hereafter simply as "abundances" un- 
less needed to avoid ambiguity) are expected to be proportional 
to the true molecular abundances (i.e., with respect to H2) given 
the stability o f CO under a number of shock conditions (e.g., 
iBergin et al.ll998) . Changes in the nature of the chemistry, how- 
ever, like from shock-processed to protostellar-wind origin, can 
potentially afifect the CO abundance and therefore the normal- 
ization factor (see Glassgold et al. 1991). 

In Table |4] we present abundance estimates with respect to 
CO for all the molecules in our survey, as derived using the val- 
ues given in Table[3] As can be seen, the highest CO-normalized 
abundances in all the outflow regimes are those of CH3OH, with 
values up to almost 10"^ in the wings. The next most abundant 
species are SO, HCN, CS, and SiO, with values of a few times 
10"^, also in the wing regimes. HC3N, para-H2CO, and HCO^, 
on the other hand, present the lowest measured abundances in 
our survey, with relative values respect to CO on the order of 
10-^ 

Given the order-of-magnitude difference in energetics and 
central source luminosity between the L1448 and 104166 out- 



flows, we start our chemistry analysis by comparing the abun- 
dance results for the two objects. This comparison is illustrated 
in Fig. Q with the LI 448/104 166 ratio of CO-normaHzed abun- 
dances for each outflow regime and for each molecular species 
detected in both sources. As it can be seen, when the abundances 
are compared in the same outflow regime, all ratios are clustered 
around 1 and show no clear deviation from that value, indepen- 
dent on the species or the outflow velocity regime. The mean 
value of the ratios in the sample is 1.2 and the rms is 0.5, which 
suggest that the abundances are very close, or equal, in the two 
sources. This is a first indication that the abundances determined 
from our survey may represent "typical" values in outflows from 
Class sources, and not be just peculiar to L1448 and 104166. 
The comparison with the LI 157 outflow abundances in section 
4.4 will further strengthen this interpretation. 

4.2. Chemical differences between the outflow regimes 

The analysis of the spectra in Sect. 3.2 already suggested that 
species like CH3OH and H2CO favor low velocities, while 
species like SiO are enhanced in the EHV gas. Table [3] shows 
quantitatively that these trends are real, and that the abundance 
of most molecules varies systematically between the different 
outflow regimes. To illustrate these abundance variations, we 
compare in Fig.[8]the CO-based abundances of several molecules 
normalized to a value of 1 in the s-wing regime. As the left 
panel shows, the abundance of CH3OH and SiO follow oppo- 
site trends across the three outflow regimes. If the abundances 
of these species are normalized to 1 in the s-wing regime, they 
differ by a factor of 20 in the EHV gas. The SO abundance, on 
the other hand, remains approximately constant across the three 
outflow regimes, indicating that this species has either constant 
abundance or that any abundance variation follows closely that 
of CO. As the figure also shows, the abundance variations be- 
tween the outflow velocity regimes are very similar in L1448 
(triangles) and 104166 (squares), illustrating again the consistent 
behavior of the chemistry in the two outflows. 

In contrast with the monotonic changes with velocity seen in 
SiO, SO, and CH3OH, the abundance of other species presents 
more complex patterns. CS and HCN, shown in the right panel of 
Fig-IU have significant (factor of a few) enhancements between 
the s-wing and f-wing regimes, but they are strongly under- 
abundant in the EHV gas. As the figure shows, HCN presents 
the most extreme change in the EHV regime, and its abundance 
drops by more than one order of magnitude respect to the f-wing 
value in both L1448 and 104166 (even assuming that the wing- 
like HCN emission we observe truly represents EHV gas). This 
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Fig. 8. Variation of CO-normalized abundances over the three outflow regimes. Values are set to unity in the s-wing regime to 
better illustrate their change across the outflow. Squares represent L1448 data and triangles represent 104166 data. Note the good 
agreement between the two outflows and the systematic differences between the species. The CS data have been multiplied by 3 to 
avoid overlap with the HCN curve. 



abundance behavior of CS and HCN indicates that the f-wing 
and the EHV regimes have rather independent chemical com- 
positions, and that molecules cannot be simply characterized as 
favoring high or low velocities across the whole outflow range. 



Further illustration of the chemical difi'erences between the 
outflow regimes is shown in Fig.|9] where we compare the LI 448 
spectra for HCN and SiO. These species follow a similar abun- 
dance increase between the s-wing and f-wing regimes, but they 
behave very dififerently in the EHV gas. To properly compare the 
spectra, we have used the SiO(2-l) line as a template to model 
the HCN(l-O) hyperfine triplet, as discussed in the Appendix 
C, and we have matched the intensities of the two lines for the 
f-wing regime. Given the similar frequencies of SiO(2-l) and 
HCN(l-O) (Table 1), their beam-dilution factors are approxi- 
mately equal, so the different velocity regimes can be compared 
directly from the spectra without additional scaling factors. As 
can be seen in the figure, the spectra match well inside the f-wing 
regime, but the two lines strongly diverge from each other in the 
EHV gas. SiO presents a bright EHV peak, while HCN shows 
a weak wing. As discussed in section 3.4, the HCN wing could 
be a continuation of the f-wing regime and not a separate EHV 
feature, so it is difficult to estimate the exact abundance contrast 
between the two species in the EHV regime. Comparing the in- 
tensities of the two spectra in the fastest second half of the EHV 
regime (to minimize wing contamination), we estimate that the 
SiO/HCN abundance ratio in the EHV gas is at least a factor 
of 20 lower than in the f-wing regime, and the factor could be 
higher if the high velocity wing in HCN is still part of the f- 
wing gas and does not belong to the EHV regime. Such a large 
difference in abundance suggests that the mechanism causing the 
enhancement of HCN in the f-wing regime is not operating in the 
EHV gas, and that a combination of several chemical processes 
is required to explain the abundances observed across the out- 
flow regimes. Before investigating possible processes behind the 
observed behavior, it is convenient to compare the abundances 
measured in the L1448 and 104166 outflows with those already 
determined for different environments. 
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Fig. 10. Abundance enhancement factors with respect to dense 
core values for the three outflow regimes of LI 448 (red squares) 
and 104166 (blue triangles). Note how all abundances, with the 
possible exception of HCO^, are significantly (order of magni- 
tude) enhanced in most outflow regimes. Note also the lower 
enhancement of species like CS and HCN in the EHV gas. 



4.3. Comparison with dense core abundances 

We first compare the outflow abundances with those expected for 
the surrounding ambient gas. Ideally, we should use the ambient 
line components in the outflow survey spectra to infer the com- 
position of the dense gas surrounding L1448-mm and 104166. 
Unfortunately, this is not possible because the ambient compo- 
nents suffer from optical depth effects and contamination from 
the reference position in the "wobbler switch" observing mode. 
As an alternative, we compare the outflow abundances with 
those of similar dense cores for which a reliable abundance de- 
termination exists. Two good sources for this comparison are the 
starless cores LI 498 and L1517B in the Taurus- Auriga molec- 
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ular cloud (to which 104166 belongs), whose composition has 
been determined in detail using multi-line radiative transfer by 
[Tafalla et al. (2004a, 2006). Like most starless cores, L1498 and 
L1517B suffer from strong molecular depletion at their centers 
due to freeze out onto dust grains. For the comparison with the 
outflow abundances, we use the undepleted abundances charac- 
teristic of the core outer parts, which most likely correspond to 
the pre-shock ambient abundances for the positions in our sur- 
vey. As the abundances in L1498 and L1517B are similar, we 
take their geometric mean and use the result as our reference 
value for the comparison. We define an enhancement factor ff„k 
for each species as the ratio between the outflow abundance and 
the reference core abundance, and present the result for each out- 
flow range in Table |4] 

Fig[TO]shows a plot of the fe„h factor for each outflow regime 
in LI 448 and 104166. As akeady mentioned, the LI 448 and 
104166 results agree with each other for all species and out- 
flow regimes for which we have measurements, and from now 
on we will ignore source to source variations and concentrate 
on a few general trends (we consider that factor-of-2 differences 
are at the limit of significance). Overall, most fe„h factors shown 
in Fig [To] are larger than unity, indicating that the outflow abun- 
dances exceed systematically the abundances in starless cores. 
As can be seen, the two wing regimes present similar patterns 
of enhancement, having the f-wing gas higher enhancement fac- 
tors than s-wing in all species but CH3OH and P-H2CO. The 
pattern of enhancements in the EHV gas, on the other hand, is 
clearly different from that in the wings, especially concerning 
CS and HCN. As discussed in Section 4.2, these species are sig- 
nificantly enhanced in both s-wing and f-wing regimes, but their 
abundance exceeds the core values only marginally in the EHV 
gas. 

As the figure shows, the enhancement factors in the out- 
flow gas span 4 orders of magnitude. HCO^ is the least en- 
hanced molecule in our survey, showing an overabundance of at 
most a factor of 2 with respect to the dense cores, and has only 
marginal variations between the different outflow regimes. A 
second group of molecules consists of P-H2CO, HC3N, CS, and 
HCN, which present abundance enhancements of about one or- 
der of magnitude in the two wing regimes. Interestingly, three of 
these species show lower enhancement factors in the EHV gas, 
while the case for HC3N is unclear since its non detection in the 
EHV regime imposes only a very weak constraint. An additional 
order of magnitude increase in the abundance enhancement fac- 
tor is found for SO and CH3OH. As seen in the previous section, 
the SO factor remains almost constant over the whole outflow 
range (at a level of about 300), while the CH3OH factor gradu- 
ally drops by about 3 from the s-wing gas to the EHV regime. 
Finally, SiO presents the most extreme enhancement factors of 
the sample. The SiO /f„/, factor increases by more than 5 be- 
tween the s-wing and the EHV regimes, and has a typical value 
in the f-wing gas of about lO'*, with the caveat that the non de- 
tection of SiO in starless cores makes this enhancement factor a 
somewhat uncertain lower limit (Ziurys et al. 1989). 

The above overabundance of most molecules in the out- 
flow gas indicates that this material has been strongly processed 
chemically with respect to the more quiescent gas around it. This 
is particularly evident in the wing component, which most likely 
arises from gas initially in the ambient core and that has been 
shock-accelerated by the outflow wind. Indeed, the pattern of 
abundances in the wing components of the outflows bears strong 
similarities with that of the LI 157 outflow, as we will now see. 
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Fig. 11. Comparison between CO-normalized outflow abun- 
dances in L1448 (red squares), 104166 (blue triangles) and 
LI 157 (open stars) for all species detected in the molecular sur- 
vey. Note the good agreement (better than a factor of two) for 
most species. See text for a discussion on the discrepancies in 
H2CO and CH3OH. For all sources, f-wing regime values have 
been used. 



4.4. Comparison with the L11 57 outflow 

The LI 157 outflow is commonly considered the reference source 
when studying the chemistry of low-mass outflows. It is pow- 
ered by a 11 Lo Class source (Umemoto et al. 1992), and 
its composition presents high molecular abundance enhance- 
ments with respect to the surround ing med ium sug gestive of a 
strong shock chemistry in action dBachill er & Pere z Gutierrez! 
1997). To compare the outflow abundances of LI 448 and 
1041 66 with those of LI 157, we have re-analyzed the data of 
Bac hiller & Perez Gutierrea(ll997h for the so-called Bl position 
(the brightest peak in the blue lobe) using the same methodology 
applied in Section 3.5 to the L1448 and 104166 data. The results 
of this new analysis are not very different from those in the origi- 
nal paper, and are presented with more detail in Appendix D. As 
the LI 157 outflow lacks an EHV component, a proper compar- 
ison with the L1448 and 104166 outflows should be done using 
the wing regime. For L1448 and 104166, we chose the f-wing 
regime, since in this range the abundance of most species re- 
mains constant with velocity, and because for molecules with 
complex spectra like HCN and CH3OH, the abundances in the f- 
wing regime are better determined than in the s-wing regime, 
which is often confused by line overlaps. For L1157-B1, we 
define an f-wing component taking the fastest half of the wing 
emission, following a similar criterion to that used to define the 
f-wing regime in LI 448 and 104166. 

Fig [TTj presents a comparison between the CO-normalized 
abundance of all the species observed both in LI 157-Bl and our 
two outflow targets. As it can be seen, the abundance of most 
species agrees better than a factor of 2 in the three outflows, be- 
ing the only exceptions P-H2CO and CH3OH, which are more 
abundant in LI 157. This general agreement between the L1448 
and 104166 abundances and those of LI 157-Bl is remarkable 
given the almost 2 orders of magnitude spanned by the values 
in the figure and by the several orders of magnitude of enhance- 
ment with respect to the ambient values that these abundances 
represent. It is also remarkable because the LI 157 outflow is 
very different in physical properties from the L1448 and 10461 1 
outflows. In addition to lacking EHV gas, LI 157 presents a com- 
bination of a low velocity range and a high CO column density 
that makes its CO column density per km s ' higher than that of 
L1448 and 104166 by factors 13 and 36, respectively. The good 
matching between abundances shown in Fig [TTj therefore indi- 
cates that outflows with different physical properties, but still 
powered by Class sources, can have very similar chemical 
compositions. It also indicates that the abundances do not de- 
pend strongly on the environment where the source has formed, 
as LI 448 and 104166 lie in large clouds (Perseus and Taurus), 
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while L1157-mm was formed in a more isolated environment. 
We take this good match as a suggestion that the abundances de- 
termined in this analysis are not peculiar to the sources in the 
sample, and that they likely represent a pattern common to other 
Class outflows. 

The similar abundances found in L1448, 104166, and 
L1157 do not imply that there is perfect chemical unifor- 
mity. L1157 itself presents factor-of-2 abundance differences 
between its several well studied p ositions (Bl versus B2, see 



s-wing 



f-wing 



EHV 



between its several well studied positions (til versus az, see 
Bachiller & Perez Gutierrez! 119971 : iRodriguez-Fernandez et al.l 



20101) . and high angular resolution obse rvations of the B 1 region 



reveal interna l composition gradients (Benedettini et alJ l2007t 
ICodella et al.l i200 9). In addition. Fig [TT] shows that the agree- 
ment between the abundances in L1448, 104166, and LI 157 is 
not perfect. Even after ignoring abundance differences smaller 
than a factor of 2 (that could arise from uncertainties in the anal- 
ysis), discrepancies remain in the abundances of P-H2CO and 
CH3OH. The factor of 5 discrepancy in the abundance of p- 
H2CO, however, should be treated with some care. For LI 157, 
only ortho H2CO has been measured, while for L1448 and 
104166, only the para form was observed, and it was neces- 
sary to assume an ortho-to-para ratio for the comparison (the 
high-temperature limit of 3 was used). Additionally, our H2CO 
abundance determination is based on a single lines, and it can 
be uncertain due to NLTE excitation in the irregularly-spaced 
H2CO energy ladder. The factor of 3 higher CH3OH abundance 
in LI 157, on the other hand, seems better established, as the 
same transitions were used in all sources. Interestingly, both 
H2CO and CH3OH are the two molecules that favor low veloci- 
ties in L1448 and 104166 (H2CO is significantly more abundant 
in the s-wing regime), and they also present a similar behavior in 
LI 157-Bl (where CH3OH is sHghtly more enhanced than H2CO 
in the slower gas). Such velocity dependence of the abundance 
could in principle explain the higher values measured towards 
LI 157, as the gas in this outflow moves slower than the gas in 
L1448 and 104166. Further observations of H2CO and CH3OH 
are needed to test this hypothesis, and a more detailed study of 
the chemistry of these two species, using a large sample of out- 
flows, will be presented in Santiago-Garcia et al. (2010, in prepa- 
ration). 

Even after considering the above differences among outflows 
(and between parts of the same outflow in the case of L 1 1 57), the 
main pattern that emerges from Fig [TT] is that of relative (factor 
of 2) homogeneity between the composition of the different out- 
flow wings. As further discussed in the next section, that com- 
position most likely arises from shock chemistry, and many of 
the observed variations are likely to result from the natural de- 
pendence of this chemistry on the shock velocity. In that sense, 
L1448, 104166, and LI 157 seem representative of a chemically 
rich early phase in the life of a bipolar flow. 

5. Nature of the outflow regimes and their chemical 
composition 

5.1. Origin of fhe observed outflow velocities 

Before discussing the origin of the different outflow regimes, it 
is worth noting that the presence of these regimes provides by 
itself useful information on the outflow velocity field. Line spec- 
tra only provide information on the radial component of the gas 
velocity field, so they do not reveal the full 3D picture of the 
gas motions. Outflow wings, in particular, can arise from a va- 
riety of velocity fields, of which there are two extreme cases 
dMarguhs & Ladalll985h : (i) the gas has a single speed and the 
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Fig. 12. Intensity ratio with respect to CO in L1448 for three 
species with different behavior across the outflow regimes. Note 
how H2CO favors the s-wing, HCN favors the f-wing, and SiO is 
mostly present in the f-wing and EHV gas (with a slight increase 
in the EHV regime). 



different radial velocities in the spectra arise solely from pro- 
jection effects, and (ii) the gas has a true distribution of speeds, 
and the observed radial velocities faithfully reflect that distribu- 
tion. The presence of velocity-dependent molecular abundances 
can help distinguish between these two extreme cases, as molec- 
ular composition, in contrast to velocity, is projection indepen- 
dent. To illustrate this idea, we present in Fig. [12] the velocity- 
dependent emission ratio with respect to CO of three species, 
H2CO, HCN, and SiO. As can be seen, each ratio depends on ve- 
locity in a different way: H2CO is more abundant in the s-wing 
regime, HCN dominates the f-wing regime, and SiO does so over 
both the f-wing and EHV gas, in agreement with the results of 
previous sections. This variety of behaviors, together with the 
close-to-constant excitation conditions derived in Sect. 3.5, in- 
dicates that the emitting gas cannot consist of a single, homo- 
geneous gas component whose distribution in radial velocities is 
only due to projection effects. The observed changes in the gas 
composition as a function of velocity imply necessarily that the 
velocity regimes we observe are intrinsically distinct, and that 
the radial velocities must represent in some way the true (3D) 
velocities of the gas. 

The idea that projection effects in a constant-velocity 
gas cannot explain the observed outflow wings also agrees 
with the high collimation of the L1448 and 104166 outflows. 
Interferometer observations of the flows show that their opening 
angle near the YSO ranges from 30° to 45°, and that th e colli- 
mation at large distances is even h igher than that (Bachille r et al.l 
' l995l:ISantiago-Garcia et al.l l2009). Projection effects from such 
a narrow range of angles seem insufficient to explain the high 
velocity span of the outflow wings (~ 50 km s"' in L1448 
and ~ 30 km s"' in 104166). Thus, in the following discus- 
sion we will assume that the observed radial velocities are pro- 
portional to the 3D gas speed of the flow. The proportionality 
factor, unfortunately, is rather uncertain, as the inclination an- 
gle of the outflows with respect to the l i ne of sight is not well 
constrained. For LI 448. 'Girart & AcordI (12001') used a 2-epoch 
measurement of the proper motions of the EHV SiO emission 
peaks to estimate an inclination angle with respect to the plane 
of the sky of 21°, which implies a true EHV velocity of about 
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180 km s '. Although possible, the projection correction im- 
plied by this measurement seems very high, as the EHV com- 
ponent in L1448 already has one of the highest radial velocities 
observed, and with the correction, L1448 becomes an extreme 
object among its class. Lacking an alternative measurement for 
104166, and to avoid a possible overcorrection factor, we will 
assume in the following discussion a representative projection 
angle of 45° for both LI 448 and 104166. We remark, however, 
that a larger correction will not invalidate the discussion below. 
In fact, it will make the problems of the chemical models dis- 
cussed to explain the EHV emission even more acute than we 
estimate. 



5.2. Two outflow wing regimes: velocity-dependent shock 
chemistry 

As discussed before, the chemical composition of the outflow 
wings in L1448 and 104166 is very similar to that of the LI 157 
outflow, which is often considered a textbook example of shock 
chemistry (Bachiller & Perez Gutierrez 1997; Benedettini et al.' 
[2007; Codellaetal. 2009; Rodriguez-Fernandez et al. 2010). 
In addition, the gas in the outflow wings is commonly in- 
terpreted as resulting from accelerated ambient gas (e.g., 
ISantiago-Garcia et al.l 12009), and an acceleration to the ob- 
served supersonic velocities necessarily implies the presence of 
a shock. Thus, it seems natural to interpret the wing abundances 
in LI 448 and 104166 in terms of shock chemistry, and to try to 
explain the composition differences between the s-wing and f- 
wing gas as a result of the velocity dependence expected for this 
type of chemistry. 

In section 4.2 (see also Fig. \T% we saw that the s-wing 
regime has lower SiO and HCN abundances than the f-wing 
gas, while it is relatively enhanced in H2CO and CH3OH. These 
latter molecules are often referred to as "primary" species, be- 
cause their abundance enhancement in outflows and warm gas 
is thought to occur via their direct relea se from the mantles 
of dust grains (e.g.. ICharnlev et aDll992h . In outflows, such a 
release is often explained as the result from the sputtering of 
the (electrically charged) dust grain s as they stream thr ough the 
neutral gas inside a C-type shock (iDraine et al.lll983h . As the 
high H2CO and CH3OH abundances in the s-wing gas indicate, 
grain mantle destruction must be occurring even at the relatively 
low velocities characteristic of this outflow regime, whose lower 
limit we estimate in 5 km s"' assuming a inclination angle of 
45° for both the L1448 and 104166 outflows. Such low veloc- 
ity starts to be uncomforta ble for grain sputtering models. The 
classical analysis by Drain eetaD (Il983h predicts that most wa- 
ter ice is returned to the gas phase due to sputtering for shocks 
faster than 25 km s ~', and although the number has been recently 
lowered slightly bv lJimenez-Serra et al.l (l2008h for the release of 
CH3OH and H2O, it is still on flie order of 20 km s ' . (A much 
lower value, around 10 km s"' was proposed for dirty ices by 
iFlower & Pineau des ForetsI (1994). but the value of the sputter- 
ing threshold used by these authors is uncertain.) Requiring that 
our observed minimum s-wing radial velocity corresponds to the 
high speed of the above models would imply that both outflows 
are less than 15° from the plane of the sky, which seems un- 
likely given the already very high (radial) velocity of their EHV 
gas. On the other hand, although it is possible that the outflow 
gas has slowed down significantly from its initial shock velocity. 
Fig. [12] shows that the SiO and HCN abundances drop in the s- 
wing regime compared to the f-wing regime. This suggests that 
only little gas from the f-wing regime has been transferred to the 
s-wing part of the flow. 



A possible solution to the observed low-velocity enhance- 
ment in CH3OH and H2O is to invoke an additional process 
of mantle disruption, like di r ect collisions between dust grain 
particles (JTielens etal .111994 . ICaselh et aO (1 19971) have shown 
that grain-grain collisions in oblique C-type shocks can lead to 
complete water ice evaporation at shock velocities as low as 
15 km s"'. These authors, in addition, suggest that the release 
of less bound species like CH3OH could occur at velocities as 
low as 10 km s"', in better agreement with our observations. 
Clearly more work is needed to understand low-velocity mantle 
disruption, which is not peculiar from the sources o f our survey, 
but se ems occur with relatively high frequency (e.g.. lGarav et al.l 
I2002h . 

In contrast with the s-wing regime, the f-wing gas is char- 
acterized by an abundance drop of H2CO (and to a less extent 
CH3OH) together with an enhancement of SiO, HCN, and CS 
(Fig. [T2I 1. This composition change occurs at a velocity of ap- 
proximately 20 km s"' for both L1448 and 104166, again as- 
suming an outflow inclination angle of 45°. As illustrated by 
the Fig. [12] the composition change between s-wing and f-wing 
regimes is relatively sudden, especially taking into account that 
after the transition, the abundance of most species remains al- 
most constant for about 20 to 30 km s"' (up to the EHV regime). 
Such a relatively sharp transition between the s-wing and f-wing 
regimes suggests that one or more additional chemical processes 
start to operate in the gas at that velocity, and that they con- 
tinue to operate for the rest of the f-wing regime. To our knowl- 
edge, no current chemical model can explain the changes we 
observe. The abundance decrease of H2CO and CH3OH, for ex- 
ample, suggest that these species are being destroyed at high 
velocities either during grain sputtering or via gas-phase reac- 
tions, while mantle sputtering and grain-grain collision mod- 
els commo nly predict a steady increase in abundance with ve- 
locity (e.g.. ICaselh et alj|1997i: fjimenez-Serra et al.ll2008l) . The 
enhancement of SiO, on the other hand, has been modeled in 
detail by a number of authors, and two alternative scenarios 
have been usually explo red: release of silicon from grain cores 
and from grain mantles (ISchilke et al.lll997t ICaselU et anil997l: 
Gusdorf et al. 2008a b). The sudden abundance increase of SiO 
near 20 km s ' could in principle be interpreted as indicating the 
crossing of the threshold for core grain disrupti on (although a 
value more 30 km s"' would have been expected. JGusdorf et al.l 
2008a). Still, the relatively constant SiO abundance in the f-wing 
gas after crossing that threshold is puzzling, since models pre- 
dict a steady increase of the SiO yield with velocity once grain 
core erosion starts. Thus, although our observed f-wing abun- 
dances are likely to result from shock chemistry, and the values 
measured towards L1448 and 104166 agree well with each other 
(and with those in LI 157), no current model seems to explain all 
the observed trends. 

5.3. The peculiar composition of the EHV component 

The composition of the EHV gas has so far been poorly known 
due to the lack of data apart from CO and SiO. The detection 
of SO, CH3OH, and H2CO, together with the less clear detec- 
tions or upper limits for HCO^, CS, and HCN, make it possible 
now to explore the chemistry of this gas component in more de- 
tail. As discussed in section 4.2, the EHV gas partly continues 
the trend of velocity-dependent abundances found in the wing 
components: SiO becomes more enhanced, SO (and probably 
HCO^) remain almost constant, and H2CO and CH3OH drop 
slightly with respect to the wing gas. CS and HCN, on the other 
hand, deviate strongly from this behavior The two species are 
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significantly (~ factor of 2) enhanced in the f-wing regime com- 
pared to the s-wing gas, but the transition to the EHV gas repre- 
sents an abundance drop of about one order of magnitude. The 
most extreme example of this behavior is HCN in L1448, due to 
its high signal-to-noise ratio (Figs.l9land[T2]). Even assuming that 
the weak, wing-like HCN feature that we observe at the high- 
est velocities corresponds to the EHV gas and is not the (more 
likely) continuation of the f-wing regime, the HCN abundance 
in the EHV component is a factor of 10-20 lower than in the f- 
wing gas. This drop is so high that any HCN enhancement over 
the ambient abundance observed in the wing regime seems to 
have disappeared in the EHV gas (FigfTOll. 

The CS and HCN abundance drop between the f-wing and 
EHV regimes represents the highest contrast seen in our survey, 
and suggests a significant chemical difference between these two 
outflow regimes. HCN and CS have in common that their en- 
hancement in shocked gas depends on the presence of carbon in 
the gas phase. HCN is mainly produced in shocks by the reac- 
tion of N with CH2, while CH2 i tself originates from the reac- 
tion of C with HJ or C+ with H2 (iPineau des Forets et al1ll990t 
[Boger & Sternberg 2005). CS, on the other hand, is enhanced in 
shocks by the reaction of SO and C, which like the previous re - 
actions, has no activation barrier dPineau des Forets et al.lll993h . 
The fact that HCN and CS can be enhanced in the wing gas but 
not in the EHV component suggests that abundance of free C in 
the EHV gas is lower than in the wing. Such a situation could oc- 
cur if most of the C in the EHV gas has been quickly locked up 
into CO, leaving behind little free C for additional reactions. If 
this were so, we would expect the EHV gas to contain a relative 
excess of oxygen. Indeed, Figs.[3]and|5]show that all molecules 
with prominent EHV feature are O-bearing species: CO, SiO, 
SO, CH3OH, and H?C O. (Note also fliat O also destroys CH2, 
the precursor of HCN, iBoger & SternberglllOOSl) . We take this 
peculiar composition of the EHV gas as evidence that this com- 
ponent is richer in oxygen (or poorer in carbon) than the other 
outflow regimes. 

The peculiar C/O abundance ratio in the EHV gas is a poten- 
tial clue to the nature of this outflow component. As discussed 
in the introduction, we can classify models of the EHV compo- 
nent as "primary wind" or "shocked ambient gas," depending on 
whether the material originates from the protostar (or innermost 
vicinity) or from the surrounding ambient cloud. Unfortunately, 
no detailed chemical models exist for the composition expected 
in either interpretation of the EHV gas. Different aspects of the 
chemistry of shocks and winds, however, have been previously 
studied, and we can use these results to explore possible dif- 
ferences between the two alternative origins of the EHV com- 
ponent. If the EHV gas represents ambient material originally 
at rest that has been shock-accelerated by a protostellar wind, 
the speed of the wind has to equal or exceed the velocity of 
the EHV component. For 104166 and L1448, the radial veloc- 
ity of the fastest EHV gas with respect to the ambient cloud is 
70 and 110 km s"', respectively, again assuming a standard cor- 
rection for an inclination angle of 45 degrees. Realistic models 
of molecular shocks in gas at densities of 10"* cm ' or higher pre- 
dict that for velocities higher than 50-70 km s ', the shock will 
dissoc iate H2 and will becom e of J (jump) type (Draine et aLl 
ri983t iLe Bourlot et all l2002h . Any shocks giving rise to the 
EHV components of LI 448 and 104166 must therefore be of J 
type, so the species observed in this outflow regime must have 
reformed from the dissociated material in the post-shock gas. 
Although models of J-type shocks predict efficient reformation 
of molecules in the postshock gas, they also predict enough UV 
radiation to photodissociate CO and keep a relatively high abun- 



da nce of neutral carbon. The detailed models of J-type shocks 
bv'Neufeld & Dal garnd(ll98 9^ show significant HCN formation 
in the postshock gas, and predict HCN column densities that are 
comparable to, and often one order of magnitude higher than 
those of SO (unfortunately, these authors do not predict CS abun- 
dances). Our observations of the EHV gas in L1448 and 104166, 
however, show a reversed situation, where the SO column den- 
sity exceeds that of HCN by one order of magnitude or more in 
the EHV component (abundances in the f-wing gas are compa- 
rable in both outflows). This overproduction of HCN (by almost 
two orders of magnitude with respect to SO) seems to pose a sig- 
nificant challenge to the interpretation of the EHV component as 
shock-accelerated ambient gas. 

To explore the alternative interpretation of the EHV com- 
ponent as a protostellar wind, we use the pioneering work of 
Glassgold et al. (1991), who studied the composition of a 150 
km s ' wind ejected qu asi spherically f rom the surface of a low- 
mass protostar. (.Panoglou et al.l l2009l have recently presented 
an alternative model of a disk wind, bu t they find no SiO pro- 
duction, in contrast with observations). iGlassgold et alJ (Il99lh 
found that although the protostellar wind starts atomic in com- 
position, it can become molecular for high enough densities 
(corresponding to mass-loss rates of a few 10' ^ Mp yr~'). Due 
to the complexity of the calculation, however, iGlassgold et al.l 
(Il99lh only considered a limited chemical network which did 
not include the chemistry of SO, HCN, or CS, so we can only 
make educated guesses on the production rates for these species. 
Given the relatively high abundance predicted for OH in the 
wind, and the lack of activati on energy for its reaction with S 
(Pineau des Forets et al.lll993 ), SO formation in the wind seems 
favorable for the same conditions for which CO and SiO are 
abundant. The possibility of HCN and CS production, on the 
other hand, is harder to estimate. It is however noticeable that 
in the denser models explored bv IGlassgold et al. ( 1991) (mass 
loss rates higher than 3 10"^ Mo yr"' and cases 2 and 3 in the au- 
thors notation), all the carbon is locked up into CO at large radii, 
and no C or C^ is available for further reactions (their Table 4). 
This lack of C suggests that HCN and CS production may not 
be favorable for the same wind models that predict abundant CO 
and SiO, in agreement with our observations of the EHV gas. 
(Note however that the model overproduces SiO, but the chem- 
istry of this species is very sensitive to the gas density and the as- 
sumed stellar temperature.) In addition to the unexplored chem- 
ical paths, the Glassgold et al. ( 1991) model contains a number 
of uncertainties that may affect the gas chemistry, like the effect 
of the higher collimation expected for a jet geometry and the role 
of the stellar UV field. Still, its chemical paths offer a promising 
mechanism to explain the observed behavior of CS and HCN in 
the EHV gas. 

More challenging for the two scenarios of the origin of the 
EHV gas is the presence of CH3OH and H2CO with abun- 
dances significantly enhanced with respect to the ambient gas 
(Sect. 4.3). Gas-pha se production of C H3OH, in particular, 
is highly inefficient (iGeppert et al.l 120061) . and any abundance 
enhancement of this species is usually interpreted as result- 
ing from the disruption of dust grain mantles, where CH3OH 
had p reviously formed via gas-grain reactions (Watanabe et al] 
1200 4). Both shocked-ambient-material and protostellar- wind 
models, however, need to rely mostly on gas-phase reactions 
to produce their molecules. In the shocked ambient-gas model, 
molecules are formed in the post-shock gas due to the disso- 
cia tive nature of t h e J-ty pe shock, while in the wind model 
of IGlassgold et al.l (|1991|) . the gas starts atomic and hot be- 
cause it originates near the protostellar surface. Despite this 
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reliance on gas phase reactions, it may be possible for either 
model to produce CH3OH and H2CO in the EHV gas by re- 
form ing it on the dust gra ins via successive hydrogenation of 
CO dWatanabe et al.ll2004 . Indeed, both the J-type shock and 
wind models predic t dust grains and abundan t atomic hydro- 
gen i n the gas phase JNeufeld & Dalga rno 1989t lGlassgold et all 
I199 I). Unfortunately, no detailed analysis of this chemical path 
has been studied so far. 

As the above discussion illustrates, much information lies 
still untapped in the chemical composition of the outflow mate- 
rial and in the variations with velocity that this composition suf- 
fers. To extract this information, detailed modeling of the chem- 
istry in both shock-accelerated ambient gas and primary winds 
from protostars and disks is needed. Also needed is a further ob- 
servational effort to characterize a large sample of outflows in 
the species found to be sensitive to gas kinematics. Chemical 
composition, due to its sensitivity to the thermal history of the 
gas and its initial conditions, offers a unique window to inves- 
tigate the so-far invisible processes that take place during the 
production and acceleration of a bipolar outflow. 

6. Conclusions 

We have carried out a molecular survey of one position in each 
of the outflows powered by L1448-mm and 104166, two Class 
YSOs in the Perseus and Taurus molecular clouds. These two 
outflows belong to a special class whose CO spectra presents a 
standard outflow wing and a discrete, extremely high velocity 
(EHV) component, and they are believed to represent one of the 
earliest phases of the outflow phenomenon. The main results of 
our study are the following. 

1. We have detected a number of molecules in the wings 
of both outflows, including CO, SiO, SO, CS, HC3N, HCO+, 
H2CO, HCN, and CH3OH. 

2. In the EHV component, we have detected definitively SO, 
CH3OH, and H2CO, as the spectra from these species present 
distinct EHV secondary peaks similar to those previously seen in 
CO and SiO. In contrast, the spectra from HCO+, CS, and HCN 
present low-level wings that overlap with the EHV regime, but 
the origin of this emission as part of a distinct EHV components 
is less clear 

3. For most molecules, we have used the intensities of multi- 
ple transitions together with a population diagram analysis to de- 
rive excitation temperatures and column densities. To compen- 
sate for differences in the telescope beam between the different 
transitions, we have estimated beam-dilution factors using high 
angular resolution CO maps of the outflows. 

4. When normalized to CO, the abundances in the L1448 and 
104166 outflows agree with each other better than a factor of 2. 
This suggests that the abundance patterns we have found may be 
common to other very young outflows. 

5. Based on chemical composition, we identify three differ- 
ent outflow regimes. Two regimes are part of the wing compo- 
nent and the third one consists of the EHV gas. The low-velocity 
wing regime presents velocity-dependent abundances for most 
species and is richer than the others in CH3OH and H2CO. The 
high-velocity wing regime has almost-constant abundance with 
velocity, and the gas is richer in species like SiO, HCN, and CS 
compared to the slower wing. Finally, the EHV gas presents a 
mixed composition: it has a very high SiO abundance but is rel- 
atively poor in HCN and CS. 

6. The outflow abundances typically exceed those of the 
Taurus dense cores by one or two orders of magnitude. On the 
other hand, the fast wing abundances are very similar to those of 



0.06 

0.04 

0.02 



0.06 

0.04 

0.02 



0.06 

0.04 

0.02 



0.06 

0.04 

0.02 



0.06 

0.04 

0.02 



0.06 

0.04 

0.02 





IRAS 04166+2706 

I (36", 72") 




(12", 36") 



■-J3licvtP-i]J-^ilf^^ 




(km s ) 



Fig. A.l. SO(32-2i) (black) and CO(2-l) (red) spectra along the 
central axis of the 104166 outflow. The CO spectra have been 
scaled down by a factor of 20, and the offsets are referred to 
the IRAS nominal position (^(72000) = 4''19'"42^.6, 5(72000) = 
-)-27°13'38"). 



LI 157, the prototypical "chemically active" flow. In agreement 
with previous work, we interpret the high abundances in the 
wing components as resulting from shock chemistry. Comparing 
the slow and fast wing regimes, we find that their composition 
differences likely originate from the sensitivity of shock chem- 
istry to the velocity of the shock, although no existing model 
seems capable of explaining all features observed in the data. 
Our observations illustrate the importance of velocity in the fi- 
nal composition of shock-accelerated gas. 

7. Concerning the EHV gas, our molecular survey shows its 
peculiar chemistry. This gas combines high abundances of O- 
rich species together with an order of magnitude drop in the 
abundance of CS and HCN with respect to the wing gas. This 
indicates that the EHV gas has a much lower C/O ratio than 
the shock-accelerated wing gas, and we suggest that it may be 
a consequence of the two outflow components having different 
physical origins. One possible explanation of this difference is 
that the EHV gas represents part of the primary wind from the 
protostar or its vicinity. Further modeling of the composition of 
this type of wind is needed to reach a definite conclusion. 
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Table A.l. EHV emission along the 104166 outflow''^ 
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Notes: (1) Velocity ranges are [-43, -26] for blue EHV and [37, 54] 
for red EHV; (2) Ofi'sets as in Fig.|AlJ (3) typical cr(I) is 0.1 K km 
S-' for CO(2-l) and 0.03 K km s"' for SO(32-2i). 



Appendix A: A survey of EHV SO emission in 
104166 

The detection of relatively strong SO emission in the EHV com- 
ponent of both LI 448 and 104166 target positions was unex- 
pected. No SO emission from EHV gas had previously been re- 
ported, so no information on its distribution and relation to other 
EHV tracers like CO and SiO was available. To better character- 
ize this SO emission, we made small maps around both sources. 
In 104166, we observed 6 positions along the outflow axis and in 
L1448 we observed three positions. This appendix presents the 
104166 data (the more limited L1448 survey is consistent with 
the 104166 results). 

Fig. lA.ll shows the SO(32-2i) and CO(2-l) spectra observed 
simultaneously along the 104166 outflow axis. As can be seen, 
EHV SO emission was detected in at least 5 of the 6 target po- 
sitions, both towards the red and blue lobes. Considering the 
limited S/N of the data and the different beam size of the two 
transitions (24" in SO(32-2i) and 11" in CO(2-l)), we find a 
good match between the SO and CO data, indicating that the 
two EHV emissions originate from the same gas. In this respect, 
we can consider the position selected for the molecular survey 
as a likely representative of the outflow as a whole. 

To further characterize the EHV SO emission in the 104166 
outflow, we study the ratio between the CO(2-l) and SO(32-2i) 
integrated intensities. Table lA. fl presents the values measured in 
our survey using the LSR velocity ranges -43 to -26 km s"' 
for the blue EHV gas and 37 to 54 km s ' for the red EHV 
gas. As can be seen, the CO(2-l)/SO(32-2i) ratio in the inner- 
most 4 positions is on the order of 20, which is consistent with 
the value derived in the molecular survey of Sect. 3. In the out- 
ermost position of each outflow lobe, however, the SO(32-2i) 
emission weakens significantly, and the CO(2-l)/SO(32-2i) ra- 
tio increases to about 40 towards the northern (blue) lobe and 
twice as much towards the southern (red) lobe. 

The close to constant CO/SO intensity ratio along the inner 
EHV jet suggests that both the SO excitation and the abundance 
relative to CO are close to constant for this outflow component 
over the central ^ 0.03 pc, and that they have values similar to 
those derived for the (8", 14") position, for which additional 
transitions of both species were observed. The increase of the 
CO/SO intensity ratio in the outer outflow, on the other hand, 
suggests that there is a drop in the SO abundance or excitation 
of the EHV gas about 0.05 pc from the outflow source. Lacking 
additional SO da ta, we use the CO and SiO observations of 
[Santiago-Garcia et alJ (l2009h to investigate this drop. The high- 
resolution CO and SiO interferometric observations showed a 
similar factor of 2-3 increase in the CO(2-l)/SiO(2-l) ratio in 



the EHV component as a function of distance from the protostar. 
The interferometric data, in addition, showed that both the CO 
and SiO emitting regions broaden with distance to the YSO, sug- 
gesting that the EHV gas is spreading laterally as it moves away 
from the pro tostar (in agreement with models of internal wo rk- 
ing surfaces. [Putrev et al]|1997l:rSantiago-Garcia et al]|2009h . In 
this case, it seems therefore likely that the excitation of SiO 
drops faster than that of CO due to its higher dipole moment, 
and that the observed change in the CO/SiO ratio is at least in 
part produced by a change in excitation. As the SO molecule has 
also a high dipole moment, it is very likely that the increase in 
the CO/SO ratio with distance from 104166 also results, at least 
in part, from a change in excitation. Multitransition observations 
are required to fully constrain the characteristics of the EHV gas 
in 104166. 

Although our observations seem to represent the first detec- 
tion of SO in EHV gas, bright SO emission in the jet-like com- 
ponents of t he HH2 1 1 and O ri-S6 outflows have been recently 
reported by iLee et alJ (1201 Ol) and IZapata et aP (l2010l) . respec- 
tively. The jet-like component of these outflows has significantly 
lower velocity than the components in L1448 and 104166, and 
their emission does not form distinct EHV features in the CO 
spectra like those shown in Fig. |2] These jet-like components, 
however, may represent a similar phenomenon to the EHV gas, 
and may owe their low radial velocity just to a location close to 
the plane of the sky. Lacking accurate measurements of outflow 
inclination angles, a comparison between the chemical proper- 
ties of these objects may help test the idea that all the jet-like 
components share a common physics (and likely origin). If so, 
chemical composition may result a more discriminant tracer of 
outflow properties than gas kinematics, and the presence of EHV 
gas in outflows may be more common among Class sources 
than currently recognized. 



Appendix B: Correction for differential beam 
dilution 

The beam size of a diffraction-limited telescope like the IRAM 
30m depends linearly on the operating wavelength, so for a given 
molecule, transitions having different rest frequencies are ob- 
served using different angular resolutions. At the CO(l-O) fre- 
quency, for example, the FWHM of the 30m main beam is about 
21 arc sec, while the ma in beam size for CO(2-l) is half that 
value dGreve et al]|1998l ). This significant change in the beam 
size with frequency, together with the small angular size of the 
regions under study, makes each transition in our survey suffer 
from a different amount of beam dilution. Such a differential ef- 
fect needs to be corrected before using the observed intensities 
to derive molecular abundances, or otherwise the high frequency 
(small beam) transitions will be over weighted, producing an 
overestimate of both the molecular excitation and the species 
column density. In this appendix we present the method we have 
used to compensate for this common observational problem. 

Properly correcting all our survey data for the effects of beam 
dilution would require making high angular resolution maps for 
each transition. In the absence of such data, we used observa- 
tions of the main outflow tracer, CO, to characterize the distri- 
bution of gas in the vicinity of our survey targets, and to es- 
timate how the observed spectrum depends on the beam size. 
To this end, we made a Nyquist-sampled map around each out- 
flow target in CO(2-l), as this transition provides an effective 
angular resolution higher than most transitions in our survey 
(a! 11"). This CO(2-l) map was used to simulate observations 
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Fig.B.l. Variation of the CO(2-l) integrated intensity as a func- 
tion of total convolving beam size (normalized to the 16" value) 
for the three outflow velocity regimes used in the molecular sur- 
vey. The intensities have been determined by convolving appro- 
priately either IRAM 30m Nyquist-sampled maps (blue squares) 
or a combined Plateau de Bure interferometer+30mdata set (red 
triangles). The dashed lines represent a set of fits with the form 
Q^g that have been used to correct the survey data for beam di- 
lution effects (see text). The error bars represent the estimated 
1-cr rms and have been calculated from the noise in the spectra 
using a baseUne fit to the Une-free channels. 



with larger beam sizes by convolving the Nyquist-sampled spec- 
tra to the required resolution. As a result, we produced a grid of 
CO(2-l) spectra with effective beam sizes ranging from about 
12" (close to the beam size for CO(2-l)) to about 30" (larger 
than the largest beam in our survey). To double check the 104166 
results (where a pointing error was suspected in the Nyquist 
map), we repeated the cal culation using the combined Plateau 
de Bure-30m data set from Santiago- Garcia et alJ (l2009h . which 
contains all the spatial frequencies, is independent of our new 
30m Nyquist-sampled map, and has an angular resolution of 
about 3". 

The results of our convolution analysis are presented in 
Fig. IB. II in the form of the CO(2-l) integrated intensities as a 
function of the total convolving beam size for the three outffow 
regimes defined in the molecular survey. As can be seen, the in- 
tensity in all regimes drops systematically with beam size, which 
is an indication that the emission is only partly resolved by the 
telescope, and that differential beam dilution truly affects our 
survey data. As can also be seen, the slope of the drop is differ- 
ent in each outflow regime, and it steepens systematically from 
the s-wing to the EHV gas. This steepening is in good agree- 
ment with the previous observation that outflow collimation in- 
creases systematically with velocity in bo th LI 448 and 104166 
(iBachiller et alJl990l:IWalla et al.l2004bl) . Also seen in Fig.iRT] 
is the good agreement between the two data sets of 104166 for 
beam sizes larger than about 14". The slight disagreement for 
smaller beam sizes partly results from the high noise in the 30m 
spectra with small beam size, as the Nyquist-sampled map is 
relatively shallow and only a few spectra contribute to the small 
beam size data points. This small beam size regime is also very 
sensitive to small pointing errors, and this seems to be the cause 
of some deviations seen in both the 104166 and LI 448 plots. 



The systematic dependence of the integrated intensity with 
observing beam size shown in Fig. lB.ll suggests that the effect of 
beam dilution can be corrected, at least to first order, using a sim- 
ple parameterization. For a symmetric Gaussian beam, the beam 
dilution factor is proportional to Qj^g if the source is point like, 
to Qj^g if the source one dimensional, and trivially equal to Qjjjg 
if the source is infinitely extended in two dimensions. It seems 
therefore reasonable that when fitting the curves of Fig. lBTl we 
use an expression of the form ^^^, where a is a free parameter 
expected to lie between and 2. This expectation is confirmed by 
the data, and as the dashed curves in Fig. lB.ll show. the analytic 
expression fits the data well inside the range of observations. 
The a coefficient in the fit is 0.6, 0.75, and 0.9 for the s-wing, 
f-wing, and EHV ranges of 104166, and 0. 1, 0.3, and 1.0 for the 
same ranges in L1448, respectively. Of particular interest is the 
close-to- 1 values derived for the EHV range in both outflows, a 
further indication that this fastest outflow regime is highly colli- 
mated. 

With the above Qj^^g fits we have converted all measured in- 
tegrated intensities into equivalent values for a 16" beam, which 
is the average size of the beams used in in our survey. The ex- 
citation and column density estimates presented in this paper 
therefore represent estimates of these parameters inside a 16"- 
diameter beam. 



Appendix C: Fitting multiple-component spectra 
with templates 

The spectra of HCN and CH3OH contain multiple components 
that partly overlap, and to calculate the integrated intensity of 
each of these components in each outflow regime, we need a 
method to disentangle them. An inspection of the data shows that 
in each multi-component spectrum, all components present sim- 
ilar shape and differ only in velocity and relative intensity, and 
this suggests that it should be possible to disentangle the outflow 
contributions by fitting the complex spectrum with a synthetic 
one made of multiple copies of the same line profile. To carry 
out such a fitting procedure, we have taken spectra from sin- 
gle lines (from CO, SO, CS, or SiO) as templates, made copies 
of them using the CLASS software, and shifted each template 
in frequency by the appropriate value for each component (as 
provided by the Cologne Database for Molecular Spectroscopy, 
see Miiller et al. 2001, 2005). By varying the relative intensity of 
each template, which is the only free parameter in the procedure, 
we have fitted the multi-component spectra, and from this fit, we 
estimate the integrated intensity of each line component in the 
three outflow regimes (only the outflow part was fitted while the 
often optically thick ambient emission was ignored). 

To identify the best fit solution, we selected a set of velocity 
ranges in the spectrum where the overlap between components 
is minor or where one component clearly dominates the emis- 
sion, and we required that the synthetic spectrum matched the 
observed integrated intensity within the noise level. For the case 
of HCN(l-O), we found that the three hyperfine components had 
relative intensities very close to their statistical weights, suggest- 
ing that the outflow emission is optically thin and that the exci- 
tation of the sublevels is close to LTE. By fixing the relative in- 
tensities of the components to their statistical weights, the only 
one free parameter left to fit the observed emission profile was a 
global scaling factor 

As shown in Fig.Hl different species present slightly different 
outflow components, reflecting their peculiar abundance pattern 
across the outflow regimes (Section 4.2). For this reason, some 
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species provide better templates than others to fit a given multi- 
component spectrum, and finding the best match is necessary a 
matter of trial and error. To test how sensitively the final inte- 
grated intensities depend on the choice of template species, we 
fitted the HCN(l-O) line from L1448 (the multiple-component 
spectrum with highest S/N in our sample) using templates from 
CS(2-1), SO(23-12), CO(2-l), and SiO(2-l). When these tem- 
plates are properly scaled to match the observed spectrum, the 
rms dispersion in the estimate of the outflow contribution to the 
s-wing, f-wing, and EHV regimes is 12%, 2%, and 14%, respec- 
tively. If these dispersions are estimates of the uncertainty in- 
troduced by the multi-component fitting procedure, we conclude 
that the added uncertainty is comparable to the uncertainty in 
the telescope calibration and in the correction for beam dilution. 
The template method, therefore, seems to provide an accurate 
way to estimate the outflow contribution in multiple-component 
profiles. 

Appendix D: Revised abundances for L1157-B1 

To properly compare the abundances of our target outflows with 
those of the L1157-B1 position, often considered as a standard 
in outflow chemis try, we ha ve re-evaluated the L1157-B1 esti- 
mates of Bachille r & Perez G utierrez ( 1997) for the species de- 
tecte d in both L1 448 and 104166. We h ave used the original data 
from lBachiller & Perez Gutierrez! ( I1997.) and we have applied to 
it the same procedure used in the analysis of LI 448 and 104166: 
correction for beam dilution, use of templates to disentangle 
overlapping components, and use of population diagram analy- 
sis f or species with more than one detect ed transition. In contrast 
with lBachiller & Perez Gutierrez! (!l997l) . who integrated the in- 
tensities across the full Une profile, we have ignored LSR ve- 
locities less than 1.8 km s"' away from the systemic velocity 
(Vlsr = 2.3 km s"') to avoid possible contamination from the 
ambient cloud. This velocity limit was set by the point at which 
the '^CO(2-l) emission drops to the noise level in the spectrum. 
As in L1448 and 104166, we have divided the line wing into two 
regimes that will be referred to as the s-wing and the f-wing (see 
limits in Table !DTb . For each outflow regime, we have derived 
beam-dependent dilution fa ctors as described in Appendix B, us - 
ing the CO(2-l) data from Bachiller & Perez Gutierrez! (!l997!) . 
and we have fitted the results with power laws having a values 
of 0.6 and 0.7 for s-wing and f-wing, respectively. 

The results of the abundance re-analysis are presented in 
Table ID. II To better determine the relatively high Te v value of 
CO ( «: 50 K), we have complemented the J=l-0 and 2-1 data of 
iBach iller & Perez Gutierrez ( 1997) with the 3-2, 4-3, and 6-5 
data from lHirano & Taniguchil (12001), also corrected for beam- 
dilution effects. For species with only one observed transition, 
we have used a default Tg^ value of 12 K as determined from 
the population diagram analysis of SiO, CS, and SO. Overall, 
the results o f our re-analysis are in good (factor of 2) agreement 
with those of !Bachiller & Perez Gutierre3 (! 19971) when the abun- 
dances are normalized to that of CO. A comparison between 
the values for the s-wing and f-wing regimes shows that, as in 
L1448 and 104166, SiO becomes relatively more abundant at 
high velocities, while CH3OH and H2CO decrease in abundance 
with velocity. The other species re-analyzed here present similar 
abundances within 20 % in the two outflow regimes. 
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Table D.l. Molecular column densities and excitation tempera- 
tures for LI 157-Bl 



MOLEC 



s-wmg' 
(cm-2) 



(1) 



T„ 
(K) 



f-wing*" 
Nr Tex 

(cm-2) (K) 



CO 


3.2 10'^ 


51 


1.1 10"^ 


52 


SiO 


1.9 10" 


12 


1.1 10" 


12 


SO 


8.1 10'^^ 


17 


2.3 10" 


15 


CS 


4.8 10'^^ 


10 


1.7 10" 


12 


CHsOH'-' 


1.0 10" 


11 


2.1 10'" 


11 


HC,N 


8.2 10" 


26 


2.5 10" 


23 


HCO+ 


5.3 10" 


12<'' 


1.7 10" 


12<3) 


o-H,CO 


6.4 10" 


12<3' 


1.6 10" 


12<3) 


HCN 


9.8 10" 


12<^' 


2.9 10" 


12® 



Notes: (1) Klsr limits are -3.5 to 0.5 km s"' for s-wing and -7.5 to 
-3.5 km s"' for f-wing; (2) The E and A-forms of CH3OH were 
analyzed separately and the column density in the table is the sum 
of the two results. The E/A ratio was found to be w 1.2 in both 
outflow regimes; (3) T„ = 12 K was assumed. 
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